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摘 要 :在 大 跨度 倒 虹 吸管 道中 有 压 阅 阀 关闭 时 通常 伴随 着 水 击 现象 , 当 水 击发 生 时 ,流体 与 国体 
的 耦合 作用 同时 存在 。 为 研究 发 生 水 击 时 ,不 同 关 阀 时间 下 的 压力 波 与 考虑 耦合 作用 时 结构 的 动 
力 响应 情况 ,以 大 跨度 倒 虹 吸管 道 与 拱 式 桥 架 结构 为 研究 对 象 , 进 行 整 体 水 力 计 算 , 基 于 流 固 耦合 
分 析 理 论 ,将 整体 水 力 模型 计算 得 到 的 水 击 压力 波 作为 流体 部 分 加 载 条 件 , 在 workbench 平台 上 进 
行 双 向 流 固 耦 合 分 析 求 解 桥架 结构 的 位 移 、 内 力 及 应 力 动力 响应 结果 。 结 果 表 明 :不 同 关 阀 时 间 工 
况 下 ,水 击 压力 波 曲线 的 变化 趋势 大 致 相同 ,结构 各 监测 点 的 位 移 时 程 曲线 同 水 击 压 力 波 曲线 规律 
基本 吻合 ;阀门 匀速 关闭 时 ,增加 闸阀 的 关闭 时 间 可 以 作为 降低 倒 虹 吸管 道 水 击 压强 的 有 效 措施 ; 
水 击 力 对 桥架 管道 的 影响 最 大 , 盖 梁 次 之 , 拱 圈 最 小 ,结构 以 发 生 顺 桥 向 动力 响应 为 主 。 
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Abstract: When the pressure gate valve is closed in the large-span inverted siphon pipeline,it is usually 
accompanied by the phenomenon of water hammer. At the same time ,the coupling effect of the fluid and 
the solid exists. In order to study the pressure wave under different valve closing time and the dynamic re- 


sponse of the structure with considering the coupling effect, the overall hydraulic calculation is carried out 
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with taking the large-span inverted siphon pipeline and the arch bridge structure as the research objects. 


Based on the fluid-structure coupling analysis theory ,the water hammer pressure wave calculated by the 


overall hydraulic model is used as the loading condition of the fluid part,and the two-way fluid-structure 


coupling analysis is carried out on the workbench platform to solve the displacement , internal force, and 


stress of the bridge structure. The results show that under different valve closing time conditions , the 


change trend of the water hammer pressure wave curve is roughly the same ,and the displacement time his- 


tory curve of each monitoring point of the structure is basically consistent with the law of the water hammer 


pressure wave curve. When the valve is closed at a constant speed ,increasing the closing time of the gate 


valve can be an effective measure to reduce the water hammer pressure; the water hammer force has the 


greatest effect on the bridge pipe ,followed by the cover beam ,and the arch ring is the smallest ,and the dy- 


namic response of the structure mainly occurs along the bridge direction. 


Key words :larsge-span inverted siphon; gate valve; two-way fluid-structure coupling; water hammer pres- 


sure wave; dynamic response analysis 


水 击 作 为 压力 管道 内 动 水 压力 变化 最 有 代表 性 
的 现象 ,是 以 一 种 压力 波 的 形式 沿 着 管 壁 传播 。 倒 
虹吸 拱 式 桥 架 上 部 倒 虹 吸 压 力 管道 ,可 能 会 由 于 出 
口 处 阀门 的 突然 关闭 ,导致 之 前 仍然 处 于 流动 状态 
的 阀门 上 游 处 水 体 突然 被 迫 停止 ,此 时 从 惯性 流动 
状态 突然 停止 的 水 体 中 蕴含 的 巨大 能 量 急需 耗 散 ， 
因此 在 水 体 中 会 产生 较 高 的 压强 (可 达 管 道 正 常 工 
作 压 力 的 几 十 倍 不 止 ) ,同时 由 于 水 体 的 可 压缩 性 又 
会 产生 沿 着 管 路 系统 来 回 传 导 的 压力 波 ,管道 系统 
内 的 压力 容易 发 生 反复 急剧 的 变化 ,往往 会 引起 管 
道 系统 的 强烈 振动 ”1 ,引发 管道 失效 ,严重 时 甚至 
会 引发 工程 事故 ,因此 ,对 阀门 关闭 过 程 中 产生 的 水 
击 现象 的 研究 显得 十 分 重要 。 

水 击 的 常用 计算 方法 有 特征 线 法 ,有 限 元 法 、 有 
限 差分 法 ,对 于 有 压 管道 水 击 作 用 的 计算 , WYLIE 
和 STREETER' 在 1974 年 提出 了 具有 开创 性 的 特 
征 线 法 ,目前 为 止 ,特征 线 法 依然 是 管道 水 击 计算 的 
主要 理论 方法 ,但 特征 线 法 是 一 维 计算 方法 ,并 未 考 
虑 管道 的 空间 效应 中 。 近 年 来 , 随 着 计算 流体 力学 
理论 (computational fluid dynamics, CFD ) 不 断 发 展 与 
计算 机 性 能 的 不 断 突 破 , 已 有 学 者 将 CFD 技术 运用 
于 大 型 流体 有 限 元 数值 仿真 分 析 中 ,发 现 数值 模拟 
结果 更 接近 实际 “。 华 临 和 GUO 等 分别 利用 
CFD 与 Fluent 软件 对 管道 中 的 阀门 关闭 进行 了 动态 
模拟 ,发 现 CFD 技术 可 以 成 功 地 用 于 水 击 现象 的 建 
模 。 李 佳 等 … 使 用 滑 移 网 格 技术 与 自 定义 函数 
(user defined functions, UDF ) 功能 分 别 对 长 直 管 道 
下 的 球阀 . 典 浆 和 闸 净 3 种 阀门 的 水 击 压 强 变化 规 


律 进行 了 模拟 研究 ,但 倒 虹 吸管 道 不 同 于 长 直 管道 。 
许 文 奇 等 “通过 建 模 分 析 研 究 不 同 关 阔 规律 对 球 
阀 流 场 和 水 击 压 强 的 影响 ,并 未 针对 闸阀 进行 详细 
研究 。 

固体 结构 在 流 场 作用 下 会 发 生变 形 或 运动 ,这 
种 变形 或 运动 反 过 来 也 会 对 流 场 产生 影响 ,导致 流 
场 发 生 改 变 , 而 流 场 的 改变 又 反 过 来 改变 作用 在 固 
体 结构 上 的 流体 荷载 的 大 小 和 分 布 ,从 而 形成 流体 
与 固体 结构 的 相互 作用 、 相 互 影 响 , 即 流 固 耦 合作 用 
(fluid structure interaction , FSI) 1。 这 种 流 固 耦 合 
作用 ,几乎 涉及 所 有 的 工业 行业 ,在 输 水 管道 中 尤为 
明显 "”。 谭 颖 汪 使 用 ADINA 软件 对 水 电站 压力 管 
道 的 水 击 进 行 了 双向 流 固 耦合 模拟 , 发现 固体 域 参 
数 对 水 击 压强 的 影响 不 可 忽视 。 于 腾 等 中 在 work- 
bench 平台 上 对 某 倒 虹 吸 工程 底部 明 钢 管 段 进行 双 
向 流 固 耦合 求解 ,结果 符合 实际 ,并 达到 工程 精度 要 
求 。 可 见 , 有 压 输 流 管道 的 流 固 耦合 问题 ,已 有 学 者 
进行 了 相关 研究 并 取得 了 一 定 进展 ,但 研究 对 象 并 
未 涉及 倒 虹 吸 桥架 结构 。 而 对 于 倒 虹 吸 结构 的 动力 
响应 方面 研究 的 相关 报道 较 少 , 陈 丹 、 王 慧 、 冯 光 伟 、 
张 炜 超 等 ”针对 倒 虹 吸 在 地 震 作 用 下 的 结构 动 
力 响应 做 了 相关 研究 ,研究 内 容 并 未 涉及 发 生 水 击 
时 的 考虑 流 固 耦 合 时 的 结构 动力 响应 。 

为 详细 研究 大 跨度 倒 虹 吸 长 直 管道 闸阀 考虑 流 
固 耦 合作 用 时 关 阀 下 的 水 击 现象 ,本 研究 借助 有 限 
元 软件 针对 大 跨度 倒 虹 吸管 道 的 闸阀 关 立 时 产生 的 
水 击 现象 展开 仿真 计算 ,获取 不 同 关 阀 时 间 的 水 
击 压力 波 的 特性 以 及 结构 的 动力 响应 。 
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1 基本 原理 


1.1 倒 虹 吸水 力 计 算 算 法 


在 进行 倒 虹 吸水 力 计算 时 ,一 般 选 用 k-e 模型 
作为 管内 水 体 流动 的 满 流 模 型 ””。 其 中 ,有 表示 
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出 3 个 基本 方程 ,分 别 为 质量 守恒 方程 ,能量 守恒 方 
程 和 动量 守恒 方程 ,可 以 用 如 下 通 式 表示 , 即 
a( Moe} PES + div(pup) = div(Tegradg) +S (7) 
式 中 : pe 量 密度 ; 9 为 通用 变量 ,可 以 表示 速度 


矢量 ww 在 x 、y 和 z 方 向 的 分 量 (w 、v 、w ) 温度 了 
等 变量 ; 厂 为 广义 扩散 系数 ; S 为 广义 源 项 ;grad 表 


和 s 表示 汕 流 耗 散 率 。A-s 汕 流 模型 可 靠 性 
高 收敛 性 好 、 内 存 需 求 低 、 应 用 面 广 , 适 用 于 无 分 
离 . 可 压缩 以 及 不 可 压缩 等 流动 问题 。 
本 研究 使 用 的 标准 -se 方程 如 下 。 


连续 方程 
Ou/ox; = 0 (1) 
动量 方程 
OU; Ou; 1 op Ou 
+W = + 一 一 十 
ot ”0%; Pp oO%; OX 
0 Ou; Ou; 
0 用 2 
人 
滑动 能 方程 (下 方程 ) 
ol ea [fe 
3 tar Pub) = 3 | pe 
(3) 
滑动 能 耗 散 率 方程 ( s 方程 ) 
ee 
ot we PU 2) | K+ 0 | 
2 
人 
Ou, 6 OU 
C = (= Ed 5 
OX,; . OX; OX, ( ) 
2 
KH: = pC 二 (6) 
式 中 : C1。、C2s 为 常 系数 ; ok 和 er。 分 别 为 上 和 s 对 


应 的 普 朗 特 数 , 均 为 常 系数 ; C, 为 经 验 常 数 。 本 研 
究 倒 虹吸 模型 参数 取 值 如 表 1 所 示 。 
表 1 ki-s 方程 模型 中 的 常数 取 值 


Tab.1 天 -sg Constant values in equation models 


C, Cie Cs Ok Ce 


3 


0.09 1.44 1.92 1.0 1.3 


1.2 流 固 耦 合算 法 基本 原理 

流 固 耦合 问题 可 由 其 耦合 方程 来 定义 ,这 组 方 
程 既 包 含 流 体 域 又 包含 固体 域 。 流 体 动力 学 在 满足 
质量 、 能 量 、 动 量 3 个 基本 守恒 定 律 的 前 提 下 ,派生 


示 梯 度 ;div 表示 散 度 , 且 
divCa) = dQ/90% + 90/0y + ga-/gz (8) 
国体 部 分 可 由 牛顿 第 二 定律 导出 , 即 
pd = Vs + 大 (9) 
式 中 : p. 为 固体 密度 ; 人 为 体积 力 矢量 ; or. 为 柯 西 应 
力 张 量 ; 4, 为 当地 加 速度 矢量 。 
流 固 耦合 必须 满足 守恒 的 基本 原则 ,表现 在 耦 
合 界 面 处 ,流体 域 与 固体 域 的 位 移 4 、 应 力 7 、 温 度 
7 热流 量 g 等 变量 也 应 相等 或 守恒 , 即 
TINe 二 Th, 


di = d。 


(10) 
df = ds 
7 = 了 
式 中 : 式 左 表 示 流 体 域 ,用 下 标 上 表示 ; 式 右 表示 固 
体 域 ,用 下 标 s 表示 。 
式 (7) ~ 式 (10) 构 成 了 流 固 耦 合 
制 方程 , 联 立 以 上 方程 组 ,并 根据 耦合 实际 情况 ,名 
定 边界 条 件 和 初始 参数 就 可 以 得 到 流 固 耦合 分 析 的 
各 项 参数 解 。 
流 国 耦合 的 解法 有 两 种 :直接 耦合 解法 和 
分 离 解 法 。 分 离 解 法 按 设 定 顺序 依次 求解 固体 或 流 
体 的 控制 方程 ,通过 流 固 交界 面 把 固体 域 和 流体 域 
的 计算 结果 相互 交换 传递 ,每 一 步 达到 收敛 后 再 进 
行 下 一 步 求解 ,直到 求 得 最 终 稳定 的 结果 。 分 离 解 
法 可 稍 作 修改 直接 利用 已 有 的 CFD 软件 ,可 用 于 求 
解 实际 复杂 的 流 固 耦 合 分 析 模 型 ,因此 本 研究 选择 


2 工程 建 模 


2.1 背景 工程 


小 鱼 坝 倒 虹吸 位 于 滨 中 引水 工程 昆明 段 ,其 总 
体 布置 图 如 图 1 所 示 , 上 接 蔡 家 村 隧洞 ,下 接 松林 隧 
洞 ,采用 拱桥 式 倒 虹 吸 跨越 赔 暑 川 。 小 鱼 坝 倒 虹吸 


了 
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水 平 长 717.4 m, 实 长 782. 2 m, 输 水 钢管 采用 3 根 
内 径 4.2 m 的 Q345C 压力 钢管 ,钢管 壁 厚 为 32 mm， 
最 大 静水 头 约 137.3 m, 倒 虹吸 设计 流量 为 3 管 总 和 
100 m’/s。, 


2.2 模型 建立 


2.2.1 倒 虹 吸 整 体 水 力 模型 建立 

本 研究 采用 Fluent 建立 倒 虹 吸 整 体 水 力 模型 ， 
根据 3 根 倒 虹 吸管 几乎 一 致 与 倒 虹 吸 结构 对 称 的 特 
性 ,选取 中 间 管 道 并 将 直径 4. 2 m 倒 虹 吸管 道 流 体 
取 中 间 面 作为 对 称 面 进行 简化 建 模 为 1/2 管 进 
行 分 析 。 

根据 进口 上 游 明 渠 底 宽 4.2 m, 水 深 7.71 m, 设 
计 单 管 流量 33. 3 m?/s, 得 到 稳 态 和 人口 边界 流速 为 
Q/4 =1.028 m/s( 0 为 设计 单 管 流量 ,4 为 明渠 截面 
积 ) ; 倒 虹 吸管 壁 均 设置 为 壁面 ,壁面 设置 为 无 滑 移 
边界 ,将 倒 虹 吸管 壁 粗糙 高 度 设 置 为 0.000 5520 ;出 
口 采用 自由 出 流 outflow 条 件 , 使 用 稳 态 压力 分 离 式 
求解 器 ,二 阶 离散 格式 ,SIMPLE 算法 ;为 保证 计算 结 
果 收 敛 性 好 ,将 收敛 残 差 标 准 均 置 设 为 10 ,设置 
完毕 最 后 初始 化 并 进行 计算 。 计 算 模 型 如 图 
1 所 示 。 


717.4 m 


门 出 口 

出 口 边 界 条 件 
为 自由 出 流 

outflow 条 件 


并 回 


桥架 入 口 


桥架 出 


1 计算 模型 
Fig.1 Computational model 

2.2.2 网 格 划 分 

采用 Fluent meshing 软件 对 1Z2 管 倒 虹 吸 流体 
结构 进行 网 格 划 分 。 为 了 提高 计算 精度 , 倒 虹 吸管 
内 水 体 几 何 模 型 均 采 用 六 面体 网 格 结构 性 网 格 , 并 
在 阀门 处 管道 连接 处 等 主要 研究 部 位 进行 网 格 加 
密 ,得 到 倒 虹 吸 整体 计算 域 网 格 数量 为 267 634 ,部 
分 网 格 划分 结果 如 图 2 所 示 。 


图 2 网 格 划分 


Fig.2 Mesh generation 
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2.2.3 瞬 态 计算 中 阀门 开 合 过 程 的 模拟 与 设置 

阀门 关闭 方式 有 变速 关 立 以 及 匀速 关 立 等 多 种 
方式 ,采用 合适 的 变速 关 阔 方法 可 以 有 效 降低 水 击 
的 最 大 压强 值 ””。 为 研究 水 击 压强 与 阀门 不 同 
关闭 时 间 之 间 的 关系 ,考虑 阀门 关闭 为 匀速 关闭 时 
的 不 利 工 况 ,设置 关 阔 工 况 为 10、20、30、60 .75s 下 
匀速 关 羡 。 

本 研究 所 分 析 的 倒 虹 吸 结构 出 口 调节 阀 为 闸 
阀 ,阀门 匀速 开 闭 运 动 由 layering 和 events 动 网 格 技 
术 同 时 实现 ,通过 自 定义 函数 (UDF) 编写 闸阀 高 度 
随时 间 变 化 的 函数 '”| 来 控制 阀门 网 格 的 运动 ,实现 
阀门 从 全 开 至 关闭 过 程 的 模拟 。 水 击 压 力 计算 一 般 
需 考虑 流体 的 压缩 性 ,因此 ,此 时 边界 条 件 改 选择 为 
压力 边界 条 件 。 
2.2.4 流 固 耦 合 分 析 模 型 建立 

对 于 底部 管 段 桥架 结构 的 流 固 耦 合 动力 响应 分 
析 ,研究 表明 有 奈 输 水 管道 水 击 振动 在 双向 流 固 耦 
合 方法 下 计算 的 结果 较为 精确 ”i 。 本 研究 倒 虹 吸 
整体 空间 尺寸 极 大 , 且 布 置 形式 复杂 ,耦合 计算 时 域 
长 ,进行 整体 计算 的 可 行 性 小 .收敛 性 差 , 且 桥架 所 
在 的 倒 虹 吸 底部 管 段 ,其 两 端 就 有 限制 管道 位 移 的 
镇 墩 , 故 本 研究 根据 倒 虹 吸 的 结构 特性 以 及 探究 桥 
架 结 构 动力 特性 ,为 主要 研究 目的 ,选取 桥架 结构 所 
在 的 底部 管 段 ,只 考虑 水 击 波 传 至 此 处 时 底部 管道 、 
桥架 和 管内 流体 的 耦合 作用 。 

本 研究 采用 线 弹 性 本 构 模 型 进行 建 模 ,固体 部 
分 模型 的 材料 属性 均 按 照 实际 工程 中 的 材料 属性 设 
置 ,桥架 结构 材料 均 采 用 C45 混凝土 ,支承 环 、 倒 虹 
吸 钢管 采用 Q345 钢材 。 模 型 中 管道 .支承 环 、 桥 架 
结构 均 采 用 SOILD186 实体 单元 ,网 格 共计 156 511 
个 单元 。 网 格 均 为 六 面体 结构 性 网 格 划 分 。 模 型 的 
边界 条 件 按照 工程 实际 情况 进行 设置 :桥架 结构 的 
排 架 底部 ` 拱 脚 采用 固 结 方式 ,约束 所 有 方向 的 自由 
度 ; 倒 虹吸 管道 在 水 流入 口 侧 和 出 口 侧 均 设 置 有 轴 
问 刚 度 30.5 kNAmm 的 伸缩 节 , 因 此 在 此 处 释放 管 
道 顺 桥 向 的 平 动 自由 度 ,以 弹 得 刚度 模拟 伸缩 节 ; 倒 
虹吸 管道 与 支承 环 之 间 采 用 Bonded 的 接触 方式 ;对 
于 支 座 刚度 均 采 用 刚性 连接 。 

管内 水 体 建 模 采 用 的 是 FLUID30 单元 , FLU- 
ID30 是 Ansys 中 的 流体 声学 单元 , 主要 用 于 模拟 流 
体 介 质 及 流 ( 声 ) 固 耦 合 问题 ,在 Workbench 中 将 管 
内 水 体 定 义 为 声学 区 域 ,其 就 会 自动 转化 为 FLU- 


投稿 网 站 :http://cjam. xjtu. edu. en 微 信 公 众 号 :应 用 力学 学 报 


470 应 用 力学 学 报 


ID30 单元 。 在 输入 管内 水 体 的 材料 属性 时 ,需要 输 
入 流体 的 材料 密度 (密度 取 为 1 000 kg/m ) 及 流体 
声速 ( 取 为 1497 m/s) ,流体 黏 性 产生 的 损耗 效应 忽 
咯 不 计 ( 黏度 Viscoty 取 0) ,流体 域 计算 模型 的 最 终 
网 格 总 数 为 128 304 个 。 本 研究 的 流 固 耦合 湿 模 态 
采用 声 固 耦 合法 计算 。 


2X10* 


六 
4X10(mm) 1 
) 


1X10: 3X10: 人 
(a) 固体 域 计 算 模型 


2X10” 4X10(mm) 


1X10’ 3Xx10" 
(b) 管内 流体 域 计算 模型 


到 3 有限 元 模型 
Fig.3 Finite element model 

为 提高 计算 的 收敛 性 及 准确 性 ,将 最 大 流速 稳 
态 运 营 时 的 水 体 计算 结果 , 即 稳 态 运营 时 的 流动 状 
态 的 水 击 波 施加 到 水 体 的 入 口 、 出 口 处 并 在 Fluent 
中 激活 管 壁 动 网 格 。 将 管道 内 壁 设置 为 流 固 耦合 交 
界面 (FSI) ,最 后 使 用 软件 自 带 的 System Coupling 数 
据 交 换 平台 将 水 体 与 桥架 结构 连接 ,workbench 中 双 


向 流 固 耦 合 设置 如 图 4 所 示 。 


图 4 双向 流 固 耦合 设置 图 


Fig.4 Two-way fluid-structure coupling setup diagram 
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3 结果 分 析 


3.1 水 击 压 强 计 算 结果 

为 了 获取 下 文 流 固 耦 合计 算 所 需 的 管道 内 压力 
波动 ,在 倒 虹 吸 计 算 域 设置 了 6 个 监测 点 Pl ~ P6， 
监测 点 布置 如 图 5 所 示 。 


5 ”监测 点 布置 图 
Fig.5 Layout of monitoring points 

经 计算 , 倒 虹 吸 模型 出 入 口 过 流 能 力 均 为 
33.2 m/s, 而 最 大 设计 流量 为 33.3 m/s, 计 算 结 果 
与 设计 流量 差 为 0.3% ,因此 计算 精度 基本 满足 要 
求 , 模 型 符合 实际 流动 规律 。 

提取 倒 虹 吸水 体 监测 点 PL ~ P6 的 压强 -时 间 曲 
线 如 图 6 ~ 10 所 示 , 将 不 同 关 阀 时 间 下 的 桥架 出 口 
监测 点 P6 的 压力 -时 间 曲 线 在 阀门 完全 关闭 的 瞬间 


10.2s = 也 
一 


0.4s 
10.4s ps 


时 间 /s 


图 6 10s 关 了 阀 压强 -时 间 图 


Fig.6 Diagram of 10 s valve closing pressure-time 


5 20.1s 
= 
5 ee 
0.3s = 
1 怀 一 一 P4 
1.0 0.3s = 


压强 /MPa 


时 间 /s 


图 7 20s 关 阔 压 强 -时 间 图 


Fig.7 Diagram of 20 s valve closing pressure-time 
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压强 /MPa 


20 25 30 3 40 45 
时 间 /s 


多 8 30s 关 阔 压 强 - 时 间 图 


Fig.8 Diagram of 30 s valve closing pressure-time 


压强 /MPa 


“40 45 50 S53 65 70 75 80 


60 
时 间 /s 
9 60s 关 了 阀 压强 -时 间 图 


Fig.9 Diagram of 60 s valve closing pressure-time 


器 


压强 /MPa 


.150255560 66970750 5 900 O05 
时 间 /s 


图 10 75 s 关 阀 压 强 -时 间 图 

Fig. 10 Diagram of 75 s valve closing pressure-time 

由 图 6 ~10 可 得 到 如 下 结果 。 

1) 不 同 关 阀 时 间 下 ,不 同 监测 点 得 到 的 压强 -时 
间 曲 线 都 有 相同 的 变化 趋势 , 即 随 着 阀门 开 度 的 减 
小 ,水 击 压强 逐渐 升 高 ,阀门 开 度 减 小 为 0 时 ,水 击 
压强 达到 峰值 ,再 逐渐 振荡 至 归 0, 且 任意 关 阀 时 间 
下 各 监测 点 压强 峰值 变化 规律 基本 一 致 。 说 明 同 一 
倒 虹 吸 结构 的 计算 结果 具有 普遍 性 ,能 够 反映 客观 
规律 。 
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2) 在 阀门 完全 关闭 前 ,各 监测 点 的 水 击 压 强 峰 
值 随 着 关 阔 时 间 增 长 而 减 小 , 呈 负 相关 ,如 监测 点 
P6 在 10s 关 了 阀 时 的 压强 值 为 2.719 MPa,75 s 关 阀 
时 P6 的 压强 值 则 为 0.404 MPa, 时 间 增 加 了 7.5 倍 ， 
而 水 击 压 强 减 小 了 6.73 倍 , 因 此 通过 增加 倒 虹 吸管 
道 阀 门 的 关闭 时 间 , 可 以 有 效 地 降低 水 击 压强 ,这 与 
谭 颖 , 许 文 奇 等 呈 的 研究 结论 一 致 ,从 侧面 验证 了 
本 研究 结果 的 合理 性 。 

3) 观察 图 11 可 得 :对 比 发 现 不 同 关 立 时 间 下 曲 
线 相位 并 不 完全 一 致 , 随 着 关 阀 时 间 的 增加 ,管内 水 
体 的 水 击 能 量 随 之 减 小 ,水 击 压力 波 曲 线 的 振荡 周 
期 缩短 ,能 量 耗 散 速率 也 随 着 加 快 。 


压强 /MPa 


时 间 步 


图 11 不 同 关 阀 工 况 下 的 相位 对 比 图 


Fig.11 Diagram of phase comparison under 


different valve closing conditions 


3.2 水 击 作 用 下 桥架 结构 的 流 固 耦 合 动力 响 
应 分 析 


管内 水 体 在 水 击 作用 下 产生 剧烈 的 压力 波动 ， 
压力 波动 传播 至 桥架 处 , 便 会 导致 管道 振动 ,再 由 此 
导致 桥架 也 产生 振动 。 由 于 倒 虹 吸 桥 架 结 构 主 体 为 
拱桥 ,两 侧 为 排 架 ,动力 响应 分 析 时 应 取 倒 虹吸 管 
道 两 侧 边 墩 盖 梁 \ 主 拱 圈 拱 脚 \1/4 跨 、 拱 顶 等 部 位 
的 关键 节点 进行 分 析 , 其 中 Cl ~ C5 为 拱 圈 监测 点 ， 
B1 ~ B8 为 排 架 盖 梁 监测 点 ,Al ~ A4 为 管道 监测 点 ， 
QD ~ 加 为 排 架 结 构 , 本 研究 的 动力 响应 均 不 考虑 静 
力 效应 , 均 使 用 水 击 荷 载 下 的 相对 值 来 描述 。 

为 方便 描述 结构 的 动力 响应 ,选取 对 比分 析 使 
用 的 监测 点 的 布置 如 图 12 所 示 。 

由 前 述 分 析 可 知 , 任 意 关 阔 时 间 工 况 下 ,水 击 压 
力 波 曲线 都 表现 出 了 大 臻 相同 的 变化 趋势 。 因 此 ， 
限于 篇 幅 , 本 研究 选取 10 s 水 击 波 对 中 间 管 道 单独 
运营 关 阀 的 工 况 进 行 流 固 看 合 分 析 , 得 到 结构 在 水 
击 下 的 动力 行为 特性 。 
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(a) 监测 点 总 布置 8 
A3 A4 
Al A2 D1~D4 
(b) 管道 监测 点 (0) 侧面 监测 点 
图 12 ”监测 点 布置 图 


Fig.12 Layout of monitoring points 


3.2.1 整体 结 

管道 和 桥架 的 位 移 云 图 和 应 力 云图 如 图 
13 所 示 。 

水 击 压 强 达 到 最 大 值 时 , 倒 虹 吸管 道 最 大 位 移 
为 7.156 mm, 最 大 等 效应 力 为 144. 78 MPa, 最 大 值 


7.156 40 max 


0.182 15 min 


2.5X10° 单位 : mm 


(a) 管道 位 移 云 图 


4.619 80 max 
4.106 50 
3.593 20 

0 


5x10" 
2.5X10’ 单位 : mm 


(c) 桥架 位 移 云图 


O 


及 | 


所 在 部 位 位 置 均 为 中 间 管 道 出 口 端 ;桥架 结构 随 着 
管道 一 起 往 蜂 中 挤 压 , 最 大 位 移 为 4. 620 mm, 最 大 
拉 应 力 为 6. 119 MPa ,最 大 压 应 力 为 5. 511 MPa ,最 
大 值 位 置 均 为 14# 蘑 梁 跨 中 ,说 明 14# 善 梁 动 力 响 应 
最 大 。 


144.780 000 max 
128.700 000 
112.610 000 


及 
32.190 000 下 
16.106000 0 5X10; 性 
ER 


0.021 752 min 


2.5X10: 单位 : mm 


(b) 管道 von-mises 应 力 云图 


6.119 00 max 


2.5X10: 


(d) 桥架 应 力 云图 


妈 13 水 击 压强 达到 最 大 值 时 结构 的 位 移 云 图 和 应 力 云图 


Fig. 13 ” Displacement nephogram and stress nephogram of the structure when the water hammer pressure reaches the maximum value 


3.2.2 位 移 结果 

图 14 为 结构 入 口 侧 和 出 口 侧 的 部 分 监测 点 时 
程 曲线 ,可 明确 说 明 此 时 结构 的 动力 行为 特征 。 

由 图 13 ~ 14 可 得 到 如 下 结果 。 

1) 图 14(a) ~ ( 缚 分 别 为 倒 虹 吸管 道 . 排 架 盖 汉 
和 桥架 主 拱 圈 入 口 侧 和 出 口 侧 监 测 点 的 位 移 时 程 
图 。 可 以 看 出 ,各 监测 点 的 位 移 时 程 曲线 波形 和 水 
击 波形 一 致 ,各 监测 点 入 口 端 和 出 口 端的 顺 桥 向 位 
移 方向 相反 ,而 横 桥 向 、 竖 向 的 位 移 方向 则 相同 。 再 


结合 


结合 图 13 ,可 以 看 出 :在 水 击 力作 用 下 ,10.4s 时 刻 ， 
水 击 波 压强 达到 最 大 正 压 ,这 时 倒 虹 吸 中 间 管 道 带 
动 着 桥架 一 起 从 两 端 往 跨 中 纵向 挤 压 ,结构 位 移 达 
到 最 大 ,然后 水 压 开始 逐渐 减 小 ,在 11.4s 达到 最 大 
负 压 ,管道 带动 着 桥架 癌 两 端 纵向 拉 伸 。 说 明 水 击 
力作 用 在 倒 虹 吸 桥架 上 时 ,首先 对 管道 产生 作用 , 随 
着 水 体 压强 的 变化 管道 在 纵 桥 向 、 横 桥 向 以 及 竖 向 
来 回 振荡 ,带动 着 桥架 一 起 来 回 振荡 ,并 随 着 时 间 增 
加 ,水 击 波 能 量 衰 减 ,结构 振幅 也 逐渐 减 小 。 
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(a) 管道 入 口 侧 控 制 点 A1、A3 处 位 移 时 程 图 


位 移 /mm 
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(c) 排 架 入 口 侧 控制 点 B1、B2 处 位 移 时 程 图 
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位 移 /mm 
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时 间 /s 
(e) 拱 圈 入 口 侧 控 制 点 C2 处 位 移 时 程 图 


图 14 
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日 
呈 
i 
要 
0 5 10 15 20 25 30 
时 间 /s 
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时 间 /s 
@ 拱 圈 出 口 侧 控 制 点 C3 处 位 移 时 程 图 


结构 人 口 侧 和 出 口 侧 的 部 分 监测 点 时 程 曲线 


Fig. 14 Time-history curves of some monitoring points on the inlet and outlet sides of the structure 


2) 根 据 时 程 图 可 以 看 出 : 倒 虹 吸 桥 架 发 生 水 击 
时 各 监测 点 的 顺 桥 向 位 移 响 应 最 大 , 竖 向 次 之 , 横 桥 
向 位 移 响 应 最 小 ,说 明 水 击 作 用 下 结构 以 顺 桥 向 位 
移 为 主 ; 出 口 端 位 移 响 应 值 大 于 入 口 端 , 且 位 移 响应 
绝对 值 B8 处 > B2 处 > B6 处 > B4 处 ,说 明 位 移 响 
应 从 桥架 两 端 往 跨 中 逐渐 减 小 ,并 且 发 生 水 击 的 点 
的 位 移 响 应 大 于 未 发 生 水 击 的 点 。 

3) 倒 虹吸 管道 位 移 值 最 大 , 排 架 盖 梁 位 移 值 次 
之 , 主 拱 圈 位 移 值 最 小 ; 主 拱 圈 1/4 处 和 拱 顶 的 位 移 
值 和 盖 梁 相 比 较 小 ,其 原因 一 是 拱 圈 刚度 较 大 ,不 易 
变形 ,二 是 水 击 波 主 要 作用 于 倒 虹 吸管 道 ,再 经 由 盖 


梁 一 排 架 一 拱 圈 的 传 力 路 线 ,所 分 担 的 力 小 于 盖 梁 。 
3.2.3 内力、 应 力 结 

为 研究 结构 的 动力 响应 ,将 桥架 结构 的 各 监测 
点 和 监测 点 所 在 截面 进行 监测 ,得 到 桥架 结构 各 监 
测 点 所 在 控制 截面 在 整个 水 击 过 程 中 的 最 大 内 力 值 
及 最 不 利 应 力 表 。 对 于 表 2 列 出 的 内 力 值 , 主 拱 圈 
监测 点 (Cl ~ C5 ) 的 面 内 指 的 是 整体 坐标 系 下 的 XZ 
平面 , 排 架 盖 梁 监测 点 (B1 ~ B8 ) 的 面 内 指 的 是 整体 
坐标 系 下 的 了 2 平面 ,而 对 于 表 3 ,管道 监测 点 (Al ~ 
A4 ) 最 不 利 应 力 指 的 是 von-mises 应 力 ,桥架 监测 点 


(Bl ~ B8 ) 最 不 利 应 力 指 的 是 第 一 主 应 力 cj( 拉 ) 和 
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第 三 主 应 力 os( 压 ) 。 
监测 点 所 在 截面 最 大 内 力 响 应 表 


Tab.2 Maximum internal force response table of the 


表 2 
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section where the monitoring point is located 


监 面 内 面 外 面 内 面 外 
测 和 前 为/ 前 为/ 名 是 。。 尊 短 / 。 毒 短 / 
大 kN ND) Nm) CN) 
Bl 180 9%5 1324 143.6 313 345.4 
B2 508.5 0.2 0.1 0.1 336.8 2218.7 
B3 25.6 274 40 22 5.3 113.8 
B4 653.1 0.2 0.2 0.1 304.1 572.0 
B5 28.7 313 352 0.0 12.2 78.0 
B6 733.4 0.2 0.0 0.1 339.8 437.9 
B7 13.6 19.2 180 171.9 314 385.2 
B8 696.7 0.4 0.1 0.0 443.4 2659.2 
C1 1264.0 30.1 0.1 0.1 563.6 12 
C2 569.9 156.2 0.1 0.1 575.9 0.8 
C3 644.2 1561 0.0 0.1 846.6 1.9 
C4 195.7 118.3 0.0 0.2 436.6 10 
C5 261.8 12.6 0.0 04 416.6 1.1 
表 3 ”监测 点 所 在 截面 最 不 利 主 应 力 响 应 表 ( 单位 : MPa) 


Tab.3 The most unfavorable principal stress response table 


of the section where the monitoring point is located( unit: MPa) 


监测 点 al(von-mises 应 力 ) 03 

0.521 

了 0. 456 

a 99.294 

过 143. 090 

前 0.662 -0.678 
> 4.769 -4.317 
全 0.223 -0.200 
0. 803 -0.785 
0.174 -0.149 
Be 1.709 -1.057 
办 0.790 一 0.802 
这 5. 826 -5.206 
全 0.081 -0.132 
2 0. 065 一 0. 108 
窗 0.069 -0.111 
全 0.032 -0.056 
Bs 0.038 -0.058 
BY 0.595 -0.495 
2 0.491 一 0.408 
和 0.558 -0.457 
本 0. 696 -0.574 


根据 表 2、 表 3 可 以 得 到 如 下 结果 。 

1) 内 力 响 应 和 位 移 响应 规律 一 致 ,水 击 作用 对 
排 架 盖 梁 内 力 影响 最 大 ,对 拱 圈 的 影响 较 小 ; 排 架 
B7 、B8 处 的 内 力 响 应 最 大 , 且 该 部 位 作为 横梁 ,其 顺 
桥 向 的 水 击 力 导 致 面 外 动力 响应 , 面 外 最 大 弯 矩 为 
2 659.2 kN . m, 面 外 最 大 剪 力 为 158 kN; 而 竖 向 的 
面 内 最 大 弯 矩 为 443. 4 kN . m, 面 内 最 大 剪 力 为 
19.2 kN , 远 小 于 其 顺 桥 向 的 响应 。 同 样 的 , 顺 桥 向 
的 力 导 致 拱 财产 生 内 力 ,因此 其 动力 响应 主要 表现 
在 顺 桥 式 向 激励 作用 下 的 轴 力 及 面 内 诲 抢 . 面 内 前 
力 。 说 明 在 水 击 力作 用 下 , 顺 桥 向 的 力 对 桥架 结构 
的 影响 更 大 。 

2 ) 最 不 利 应力 响 应 同样 和 位 移 响应 规律 一 致 ， 
均 体现 出 了 上 述 水 击 作 用 下 结构 的 几 个 特点 :出 口 
侧 比 入 口 侧 响应 大 、 桥 架 两 端 比 蜂 中 响应 大 、 排 架 盖 
梁 比 主 拱 圈 响 应 大 ,因此 主 应 力 结果 也 能 非常 直观 
的 体现 出 各 个 监测 点 的 受 力 情况 。 


4 结 论 


本 研究 以 流 固 耦合 分 析 理 论 为 基础 ,将 整体 水 
力 模 型 计算 得 到 的 水 击 压力 波 作为 流体 部 分 加 载 条 
件 , 以 10s 关 阀 时 间 工 况 为 例 , 在 workbench 平台 上 
进行 双向 流 固 耦 合 分 析 结 构 的 响应 情况 ,得 出 如 下 
结论 。 

1) 不 同 的 关 阀 时 间 对 管道 的 水 击 规律 没有 影 
响 , 即 任意 关 阀 时间 工 况 下 ,水 击 压力 波 曲线 的 变化 
趋势 大 致 相同 。 并 且 在 流 固 耦合 工 况 下 ,桥架 结构 
各 监测 点 的 整个 水 击 过 程 的 位 移 趋 势 基本 和 水 击 波 
曲线 规律 吻合 。 

2) 阀门 关闭 时 间 越 短 ,使 管道 内 的 压力 变化 越 
剧烈 ,管道 内 各 监测 点 产生 的 水 击 压强 越 大 , 越 容 易 
对 管道 造成 破坏 ;而 随 着 关 阀 时 间 的 增加 ,管内 水 体 
的 水 击 能 量 随 之 减 小 ,水 击 压力 波 曲 线 的 振荡 周期 
缩短 ,水 击 压力 波 的 传播 过 程 中 的 能 量 耗 损 速率 加 
快 ,因此 增加 倒 虹 吸管 道 闸 阀 的 关闭 时 间 可 以 作为 
降低 倒 虹 吸管 道 水 击 压强 与 桥架 结构 在 水 击 作用 下 
的 动力 响应 的 有 效 措施 。 

3) 水 击 波 为 纵 桥 向 动 荷 载 ,对 结构 的 影响 也 主 
要 体现 在 顺 桥 向 的 动力 响应 上 ,具体 表现 为 :水 击 作 
用 下 ,结构 以 顺 桥 向 位 移 为 主 , 竖 向 次 之 , 横 桥 向 最 
小 ;在 水 击 力 的 激励 下 , 排 架 盖 梁 与 主 拱 圈 顺 桥 向 的 
内 力 响 应 明显 大 于 其 他 方向 ,因此 在 设计 时 要 着 重 
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考虑 桥架 结构 的 顺 桥 向 刚度 。 

4) 发 生 水 击 时 ,管道 由 于 动 水 压力 的 作用 ,会 在 
XX 轴 方向 上 来 回 振荡 ,振荡 幅度 随 着 时 间 推 移 逐 渐 
变 小 ,而 支承 管道 的 排 架 结 构 , 由 于 刚度 不 足 不 能 
效 限制 管道 的 变形 , 便 随 着 管道 一 起 来 回 摆动 , 当 摆 
动 传导 致 主 拱 圈 时 ,水 击 能 量 已 被 消耗 大 半 , 因此， 
管道 动力 啊 应 > 排 架 盖 梁 动力 啊 应 > 主 拱 圈 动力 响 
应 。 在 结构 设计 时 要 着 重 考虑 桥架 结构 的 管道 的 受 
力 性 能 ,在 管道 天 阀 时 要 时 刻 监控 结构 管道 的 受 力 
情况 。 

5 ) 在 倒 虹 吸 结构 实际 运营 过 程 中 ,可 能 会 出 现 
单 管 运营 的 工 况 。 而 在 单 管 水 击 作用 下 ,桥架 结构 
出 口 侧 与 管道 连接 的 盖 梁 构件 会 发 生 大 幅度 的 不 均 
匀 变 形 , 受 力 极为 不 利 ,相应 的 排 架 底部 也 会 产生 较 
大 拉 应 力 , 建 议 这 些 部 位 在 设计 时 也 应 着 重 考虑 ,并 
考虑 可 以 在 倒 虹吸 管道 支承 与 桥架 盖 梁 之 间 布 置 减 
隔 震 支 座 , 改 变 顺 桥 向 刚度 ,使 其 在 水 击 力作 用 下 吸 
收 缓冲 大 量 外 荷载 能 量 , 将 其 作为 有 效 阻隔 上 下 部 
结构 振动 的 措施 , 以 降低 桥架 结构 在 水 击 力 下 的 动 
力 响 应 , 支 座 具体 刚度 值 应 结合 结构 的 其 余 动力 特 
性 及 实际 情况 选取 。 

6) 在 水 击 力作 用 下 ,水 击 波 为 纵 桥 向 动 荷 载 , 结 
构 的 表现 形式 为 来 回 振 动 , 因 此 后 续 应 进一步 考虑 
分 析 可 能 发 生 的 水 击 共振 现象 及 水 击 疲劳 现象 ; 

7) 本 人 研究 均 为 有 限 元 仿真 分 析 , 结 果 均 是 基于 
理论 分 析 所 得 到 的 ,并 无 实际 试验 数据 支撑 ,下 一 步 
计划 将 双 管 运营 工 况 .3 管 运营 工 况 以 及 水 的 流速 、 
关 立 方案 与 天 阀 规律 等 其 他 应 考虑 的 影响 水 击 压 强 
大 小 的 因素 考虑 进 有 限 元 模型 中 ,并 结合 相关 规范 
等 具体 数值 进行 分 析 比 较 , 以 得 到 更 具有 代表 性 、 适 
用 性 的 结果 。 
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